Heavy Element Nucleosynthesis 


A summary of the nucleosynthesis of light elements is as follows 


^He Hydrogen burning 

3He Incomplete PP chain (H burning) 
7H, Li, Be, B Non-thermal processes (spallation) 
HN eGo Ne O CNO processing 

HO Helium burning 

180, ?2Ne a captures on 4N (He burning) 
20Ne, Na, Mg, Al, ??Si Partly from carbon burning 

Mg, AL Si, P, S Partly from oxygen burning 


Ar, Ca, Ti, Cr, Fe, Ni Partly from silicon burning 





Isotopes heavier than iron (as well as some intermediate weight iso- 
topes) are made through neutron captures. Recall that the prob- 
ability for a non-resonant reaction contained two components: an 
exponential reflective of the quantum tunneling needed to overcome 
electrostatic repulsion, and an inverse energy dependence arising 
from the de Broglie wavelength of the particles. For neutron cap- 
tures, there is no electrostatic repulsion, and, in complex nuclei, 
virtually all particle encounters involve resonances. As a result, 
neutron capture cross-sections are large, and are very nearly inde- 
pendent of energy. 





To appreciate how heavy elements can be built up, we must first 
consider the lifetime of an isotope against neutron capture. If the 
cross-section for neutron capture is independent of energy, then the 
lifetime of the species will be 
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For a typical neutron cross-section of (o) ~ 10 7? cm? and a tem- 
perature of 5 x 108 K, c, ~ 10? /N, years. 


Next consider the stability of a neutron rich isotope. If the ratio of 
of neutrons to protons in an atomic nucleus becomes too large, the 
nucleus becomes unstable to beta-decay, and a neutron is changed 
into a proton via 


(Z, A+1) — (Z+1,A+1) +e° + % (27.1) 


The timescale for this decay is typically on the order of hours, or 
~ 10^? years (with a factor of ~ 10° scatter). 





Now consider an environment where the neutron density is low 
(N,, ~ 10? cm~?), so that the lifetime of isotope A, Z against neu- 
tron capture is long ~ 10* years. This happens during the helium 
thermal pulses of an AGB star, where, due to mixing of shell mate- 
rial, the reactions 


136 + 4He + 160 4; 


14N + tHe — 150 + He —> ??Ne + He — ?Mg +n 





release free neutrons at a slow rate. In time, the isotope will capture 
a neutron, undergoing the reaction 


(Z, A) +n — (Z, A1) +7 


If isotope Z, A+1 is stable, then the process will repeat, and after a 
while, isotope Z, A+2 will be created. However, if isotope Z, A+1 is 
unstable to beta-decay, the decay will occur before the next neutron 
capture and the result will be isotope Z+1, A+1. This is the s- 
process; it occurs when the timescale for neutron capture is longer 
than the timescale for beta-decay. 








Next, consider a different environment, where the neutron density 
is high, N, ~ 107? cm~’. (This occurs during a supernova explo- 
sion.) At these neutron densities, the timescale for neutron capture 
is of the order of a millisecond, and isotope Z, A+1 will become 
Z,A+2 via a neutron capture before it can beta-decay. Isotope 
Z, A+2 (or its daughter isotope if it is radioactive) would therefore 
be an r-process element; its formation requires the rapid capture of 
a neutron. 





Whether an isotope is an s-process or r-process element (or both) 
depends on both its properties and the properties of the isotopes 
surrounding it. In the example above, suppose the radioactive iso- 
tope Z—3, A+1 is formed during a phase of rapid neutron bom- 
bardment, and suppose that its daughters Z—2, A+1, Z—1, A+1, 
and Z, A+1 are all radioactive as well. Isotope Z-F1,A--1 would 
then be the end product of the r-process, as well as the product of 
s-process capture and beta-decay from Z, A. On the other hand, if 
isotope Z—1, A+1 is stable, then the beta-decays from Z—3, A+1 
wil stop at Z—1, A+1, and element Z+1,A+1 would effectively 
be “shielded” from the r-process. In this case, element Z+1, A+1 
would be an s-process isotope only. Isotope Z—1, A4-1 would, of 
course, be an r-process element, and, if isotope Z—2, A were unsta- 
ble, it would be “shielded” from the s-process. 





A third process for heavy element formation is the p-process. These 
isotopes are proton-rich and cannot be formed via neutron capture 
on any timescale. P-process isotopes exist to the left of the “valley 
of beta stability" and thus have no radioactive parent. These iso- 
topes are rare, since their formation requires overcoming a coulomb 
barrier. Their formation probably occurs in an environment similar 
to that of the r-process, i.e., where the proton density is extremely 
high, so that proton captures can occur faster than positron decays. 
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(Neutron-rich matter) 
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A characterization of a portion of the chart of nuclides showing 
nuclei to the classes of s, r, and p. The s-process path of (n, y) 
reactions followed by quick /-decays enters at the lower left and 
passes through each nucleus designated by the letter s. Neutron- 
rich matter undergoes a chain of 6-decays terminating at the most 
neutron-rich (but stable) isobars; these nuclei are designated by the 
letter r. The s-process nuclei that are shielded from the r-process 
are labeled s only. The rare proton-rich nuclei, which are bypassed 
by both neutron processes are designated by the letter p. 
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The valley of stability. 


Isotopic Ratios 


The abundance of every isotope A, Z can therefore be written 
N(Z, A) = N,(Z, A) + N,(Z, A) + Np(Z, A) e NS(Z, A) + N,(Z, A) 


By knowing the neutron capture cross-sections of the isotopes, and 
by knowing which isotopes are shielded, (i.e., which are solely s- 
process and which are solely r-process) the pattern of elemental 
abundances can be used to trace the history of all the species. 








Clearly, the relative abundance of a given s-process element depends 
on its neutron capture cross-section. If the capture cross-section is 
low, then nuclei from the previous species will tend to “pile-up” 
at the element, and abundance of the isotope will be high. Con- 
versely, if the isotope has a large neutron-capture cross-section, it 
will quickly be destroyed and will have a small abundance. S-process 
abundances can therefore be derived from differential equations of 
the form 





dN 
dt 

The starting point for these equations is °°Ni (although neutron 
captures actually start during helium burning, since the reactions 
13C(a,n)'®O, O(a, n) Ne, Nefa, n) Mg and ??Ne(a, n)??Mg 
all create free neutrons). The endpoint of the s-process is ?99pi, 
since the next element in the series, ^! Bi a-decays back to ?99Pb. 
Thus, the entire network looks like 


= —GANA t GA-1NA-1 

















dNs6 
— A 
di O56V56 
dN 
a ——GANXd-Na 57 < A < 209; A Æ 206 
dN. 
299 eae Node e aa N ao (27.2) 


dt 


with the boundary conditions 


Nse(t =0) A=56 


Na(t—0) = 
al ) m A> 56 


These equations can be simplified locally by noting that the neutron 
capture cross-sections for adjacent non-closed-shell nuclei are large 
and of the same order of magnitude, so that 


CANA ~ o4-1N 4-1 (27.3) 


This is called the local approximation; it is only good for some ad- 
jacent nuclei. Over a length of A ~ 100, abundances can drop by a 
factor of ~ 10. 


The abundance of r-process isotopic ratios in nature can be esti- 
mated by differencing the observed cosmic abundances from that 
predicted from the s-process. Part of the r-process isotopic ratios 
may be due to processes associated with nuclear statistical equilib- 
rium under extremely high densities (p > 10'° gm-cm ?. (Under 
these conditions, the Fermi exclusion will force electron captures in 
the nucleus, and cause a build up a neutrons.) However, r-process 
isotopes can also be built up in more moderate densities, if the 
neutron density is high (~ 107? cm~?). Under these conditions, a 
nucleus will continue to absorb neutrons until it can capture no more 
(i.e., until another neutron would exceed the binding energy). At 
this point, the nucleus must wait until a beta-decay occurs. Thus, 
the abundance of an isotope with charge Z will be governed by the 
beta-decay rates, t.e., 





dNz 

EUN — Az—1Nz_1(t) m Az Nz(t) 
where Az is the beta-decay rate at the waiting point for charge 
Z. (In practice, under these circumstances, one must also take into 


account proton and a-particle captures as well, since these reactions 
will also add charge to neutron-rich material. However, because of 
the electrostatic repulsion, these capture cross-sections will be much 
lower.) R-process nucleosynthesis stops at Z — 94 (N — 175) where 
nuclear fission occurs to return two additional seeds to the capture 
chain. 
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Measured and estimated neutron-capture cross sections of nuclei on 
the s-process path. ‘The neutron energy is near 25 keV. The cross 
sections show a strong odd-even effect, reflecting average level den- 
sities in the compound nucleus. Even more obvious is the strong 
influence of closed nuclear shells, or magic numbers, which are asso- 
ciated with precipitous drops in the cross section. Nucleosynthesis 
of the s-process is dominated by the small cross sections of the 
neutron-magic nuclei. 
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The solar system o N, curve, i.e., the product of the neutron capture 
cross sections for kT' — 30 keV times the nuclide abundance per 
10° silicon atoms. The solid curve is the calculated result of an 
exponential distribution of neutron exposures. 


10 Long-time solution 
T, = 1.0 


n, 10?* 


Cycle time — 4.9 sec 
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Calculated r-process abundances after long times. The nuclei have 
been driven around a fission cycle until the abundance distribution 
becomes characteristic of a steady state. The abundance of each 
nucleus grows exponentially with an efolding time of 4.9 sec. This 
time, calculated on the basis of -decays rates, depends upon the 
neutron density and the the temperature. 
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The abundances of the elements in the solar system. The dots 
represent values obtained from the strengths of solar absorption 
lines, while the line is based on chemical evidence from the earth 
and meteorites. 


Nucleocosmochronology 





It is possible to obtain an independent estimate of the age of the 
Solar System (and, indeed the universe) from the study of radioac- 
tive species (and sometimes their daughter elements). The idea is 
simple. Consider an radioactive species, ? present in an object that 
was born at time T = t in the universe. Prior to time t, this element 
was created in various events, and had a production rate P;p(t). 
The abundance of that species today (at time t = to) will be 





t 
N = Bron p(t) eM) dr 
0 


which simplifies to 
t 
N = Beto | era: (27.4) 
0 


where A is the decay rate of the species. Obviously, calculations of 
the production rate can be treacherous. However, if one takes the 
ratio of two elements which are presumed to be created together 
(say, in r-process events), then most of the unknowns cancel out. 


So 





N; F, e ^ito in p(T) eit dr 


N, | BP,e-^sto fi p(T) eT dr 





Moreover, if the second element is stable (with A, = 0), then 


Ni = Fi dito Jo p) eT dr 
Ns Hh F p(t) dr 
Ta 
x = (1 — A; (to — (t))) (27.5) 


In other words, by observing the ratio of two species, and knowing 
the ratio of their production ratio in a star, one can obtain to — (t), 
the mean age of the elements. 


Nucleus Half-Life (yrs) 
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Long-lived Galactic Chronometers 
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Daughter Decay Mode 
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Process 


Si burning 
Si burning 
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Obviously to do this experiment, one needs to work with elements 
that are formed in the same event. In the solar system, one can re- 
strict the experiment to r-process only isotopes. However, it is not 
possible to separate out the different isotopes via stellar absorp- 
tion lines. (This is not strictly true for the lightest species such as 
deuterium and CNO, but for heavier elements it is certainly true.) 











However, elements heavier than ?'9Bi can only be made with the 
r-process (since *!°Bi a-decays back to 7°°Pb). 


The two most popular elements to perform this experiment with 
is Uranium and ???Th (half-life 13.9 Gyr) and ?°°U (half-life 4.51 
Gyr). But these lines are weak!) For example, here are the results 
from the extreme metal-poor halo star BD--17 3248 (from Cowan 
et al. 2002): 











Chronometric Age Estimates for BD+17° 3248 


Element Age (Gyr) 
Pair Predicted Observed Solar Best Limit 


Th/Eu 0.507 0.309 0.4614 10.0 > 8.2 
Th/Ir 0.0909 0.03113 0.0646 21.7 >14.8 
Th/Pt 0.0234 0.0141 0.0323 10.3 > 16.8 
Th/U 1.805 7.413 2.32 >13.4 >11.0 
U/Ir 0.05036 0.0045 0.03869 >15.5 > 13.5 
U/Pt 0.013 0.0019 0.01846 >12.4 > 14.6 




















a F BD+17°3248 
—0.8r 

ECCL ---- HD 122563 
o [s f f i l i à i i l 
A 10 r T r r T r r PR T 
RE 

c 

© 

jaw 














0.9 f ES \ | 
[ | | 
eae V e.e... BD+17°3248 ] 
OT Th Il ] 
C 1 L L | - i : E 14 
4019 4020 


Wavelength (A) 


Fic. 8.—Observed Keck/HIRES and synthetic spectra in region 
surrounding the Th n 4019.12 line, in the manner of Fig. 3. In the bot- 
tom panel, the synthetic spectra were computed for abundances 
loge(Th) oo, —1.58, —1.18, and —0.78. 
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Fic. 9.—Spectra of BD +17°3248 in region surrounding the U n 
3859.57 line. In both panels the observed Keck/HIRES spectrum is com- 
pared to four synthetic spectra computed for abundances log e(U) = —oo, 
—2.60, —2.10, and —1.60. (Note that some of the lines overlap.) The only 
difference between the two panels is the vertical scale, which is expanded in 
the bottom panel for easier assessment of the extremely weak U r feature. 


